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Abstract 

We consider directed percolation with an absorbing boundary in 1 + 1 and 2 + 1 
dimensions. The distribution of cluster lifetimes and sizes depend on the boundary. 
The new scaling exponents can be related to the exponents characterizing standard 
directed percolation in 1 + 1 dimension. In addition, we investigate the backbone 
cluster and red bonds, and calculate the distribution of living sites along the ab- 
sorbing boundary. 
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1 Introduction 



Directed percolation (DP) is the common name for spreading processes with 
active ("live") sites and an absorbing ("dead") state. The process of directed 
percolation is more precisely defined on a lattice where at each time step 
a site can become alive with probability p, if and only if at least one of its 
neighbors was alive at the previous time step [1]. For p the process rapidly 
terminates whereas for p ^ 1 the live sites spread without limit. There exists 
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a critical value p = Pc such that the distribution of clusters originating from 
an initial live site is a power law. The probability that there are living sites 
at time t after the introduction of one living site in a infinite lattice of dead 
sites at time t = is [2] 

P{t)^t-', 8^(5/v\^. (1) 

Here j3 is the order parameter exponent which determines the probability 
{e = p — Pc) for the cluster to survive to infinity, and i/y and i'± are the 
correlation length exponents along (^y ~ e'"^^) and perpendicular ~ e"''-^) 
to the time direction of the DP process. Both /3 and z/y enter in Eq. (1) which 
reflect the fact that the final survival of the cluster is determined when the 
cluster size reaches the correlation length ; the probability to reach ^y scales 

as ~ Cy '^^'^'^ in accordance with Eq. (1). 

Directed percolation is known to be equivalent to Reggeon held theory [3,4] 
which describes the evolution of a density p of five sites for a large class of pro- 
cesses. An understanding of facets of directed percolation have implications for 
the understanding of a number of widely different problems in physics, includ- 
ing the dynamics of chemical reactions and catalyzers [1,2], contact processes 
[5], spatio-temporal intermittency [6-8], self-organized criticality [9-12], and 
directed polymers [13]. 

In the present paper we consider directed percolation with an absorbing bound- 
ary [14,15]. In the context of spatio-temporal intermittency [6-8], the process 
of directed perclation can be viewed as lateral convection of intermittent spots 
along a boundary which is absorbing (i.e., the boundary layer has laminar 
ffow). Thus the active sites represent intermittent spots that originate from 
one point on the boundary due to a impurity or an edge. The deacy of the 
turbulent spots into the absorbing (laminar flow) could possibly be connected 
to directed percolation with a wall. Another more traditional picture is a cat- 
alytic process that is initiated at the edge of the catalyzer. As a result, the 
propagation of active live sites can then only propagate into the system. We 
consider very slow activation on the boundary: the activity in the system is 
supposed to be entirely determined by the propagation of one initial spark at 
the boundary at time zero. 

We consider time-directed percolation in a geometry as shown in figure 1, and 
will refer to the d transverse direction as space {x) and to the longitudinal 
direction as time (t) . The boundary absorbs the parts of the DP cluster which 
in an open geometry would have continued spreading on the other side. We 
will first discuss the 1 + 1-dimensional {d = 1) directed site DP on a square 
lattice where pc = 0.705485, /S = 0.2765, z/y = 1.7338, and = 1.0968 [15- 
18]; in the following, the subscript '1' denotes critical exponents for DP with 
an absorbing boundary. Recently, series expansions for 1 -|- 1-dimensional bond 
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Fig. 1. a) Schematic drawing of a DP cluster, b) Schematic drawing of the same DP 
cluster but now with the absorbing boundary present. One notices the existence of 
a different distribution of voids compared to the case (a). The time axis is directed 
downwards. Note that in our implementation of the absorbing boundary wc allow 
for activity on the boundary; not allowing for this simply amounts to shifting the 
absorbing boundary one lattice spacing to the left. 



2 Critical behavior 



First we notice that there remain a critical value pc for DP with an absorb- 
ing boundary, and that this value is independent of the presence of such a 
boundary. In figure 2 we display the meandering (squared) x'^{t) and the mass 
m{t) of the cluster measured for a value of p close to Pc = —0.0015). 
The meandering scales as ~ t^^^/'^w^ with t determined by ^y. The mass 
scales as ^||^±e^ ~ gj^^j both these quantities scale in the 

same way as if there was no boundary. Thus, both the meandering exponent 
X = = 0.6327 and the mass scaling exponent 1 + (z/^ — l3)/v\\ = 1.4732 

are independent of the absorbing boundary. We also measure the lifetime 
distribution and average lifetime (see below) and conclude that the three ex- 
ponents (/?, and u±) for DP are still present for DP with the presence of 
a wall. This is consistent with the fact that DP clusters with a wall consist 
of subsets of DP clusters without a wall. Therefore the correlation length and 
thus the correlation length exponents should be identical in the two cases. 
Further, for the subset of clusters that survive to a given time t, there exist a 
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Fig. 2. Cluster mass, m{t), and meandering, x^(t), as function of time in 1 + 1 
dimensions. We obtain the mass scaling exponent 1 + (z^_l — l3)/u\\ = 1.475 ± 0.005, 
and the meandering exponent 2% = 2v_\_/v\^\ = 1.265 ± 0.005, in agreement with the 
values for DP without an absorbing boundary. 



connected leftmost path closest to the wall that merge the starting point with 
a point that survive to time t. To the right of this path the cluster is identical 
to a normal DP cluster with the same scaling behavior of the density. The 
fact that the scaling of the mass of a DP cluster with and without a wall are 
identical show that for the subset of clusters that survive to a given time, the 
active points within these clusters are on average so far from the wall that the 
density is not influenced by the presence of the wall. 

Next, we investigate the distribution of lifetimes in the case of an absorbing 
boundary. In figure 3 we show the probability that there are still living sites 
at time t. We find that the presence of an absorbing boundary changes the 
exponent of the lifetime distribution from the one given by (1) (i.e., 6 = 0.1594) 
to 

Pi{t) - t-^\ 5i = 0.420 ± 0.005. (2) 



We have 6i = Pi/u^, where the order parameter exponent f3i determines the 
probability e'^^ for the cluster to survive to infinity. This yields the value 
Pi — 0.728 ± 0.008, and we conclude that /3i ^ /3. Note, as mentioned, that P 
is still needed for DP with a boundary in order to describe the scahng of the 
cluster mass near p — Pc- From the full scaling behavior 

Pi(i)=t-^V(t6^"), (3) 
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we obtain the average lifetime 

{t)^e-^\ ri = i/||(l-5i) = i/||-/?i. 



(4) 



In Ref. [15] they obtain ti = 1.000... and conjecture that Ti is exactly equal 
to unity. This leads to Pi = u\\ — 1, and 6i = 1 — l/u^ = 0.4232 in agreement 
with our result. We also measure (t) directly (Fig. 4) and obtain the estimate 
Ti = 1.000 ± 0.005, yielding an estimate for f3i in agreement with the one 
above. Even though we do not know of an analytical argument for ri = 1, we 
can obtain an upper bound by the following heuristic argument: if we assume 
that the statistics of the DP cluster are described by the statistics of the "outer 
envelope" which scales as we obtain that the first return of the envelope 
back to the absorbing boundary determines the cluster lifetime distribution. 
The probability that the envelope returns to "zero" (i.e., to the absorbing 
boundary) is accordingly given by the scaling of the first return of a fractional 
Brownian motion with exponent 

This yields the survival exponent 5' = (2 — x) — 1 = 1 — and thus the 

estimate r[ — u±, only a httle larger than the correct value. The reason that 
the 'envelope' argument is only approximate is due to the fact that clusters 
can spontaneously die out before the envelope returns to the boundary. 

Our results for the probability Pi (s) to have a cluster of size larger than s show 
that it follows a power law with an exponent — 1 = 0.283 ±0.005. This value 
can easily be explained in terms of the lifetime distribution investigated above. 
A cluster of lifetime t will have a size (mass) s of the order s ~ f^+i<i'^±-l3)/'^\\ ^ 
where the density exponent /5 enters (and not the order parameter exponent 
/3i associated with a seed on the boundary). Transforming the probability 
distribution (3) for Pi{t), we thus obtain 

Pi{s) = s-(--^)^(5e^ii+'^^--0, T, = 1 + —-^ -. (6) 

This yields = 1.287 in nice agreement with our numerical result in 1 + 1 
dimension. Compared to the value 1 + /3/(i/|| + iy± — P) — 1.108 for DP with no 
boundary we find as expected that the boundary suppresses the cluster sizes. 
Prom Eq. (6) we find the average cluster size 

(s)~e-^\ ^,^jy^^+du^-p-Pi, (7) 
which replaces the usual hyperscaling relation when there is no wall present. 
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Fig. 3. Probability Pi{t) that there are still live sites at time t in 1 + 1 dimensions. 
We obtain a power law with exponent 6i = 0.42 it 0.01. We also display the distri- 
bution Pfirst(i) of first return along the absorbing boundary and obtain a power-law 
distributions with exponent rgrst = 1-577 it 0.005. 

In standard DP the cluster mass and the order parameter scale with the same 
exponent and there are three independent exponents. In DP with an absorbing 
boundary, we find that the mass and the order parameter scale with different 
exponents. One of these exponents can however be related to the other 
exponents by assuming that Ti — 1 and we are left with the same three 
independent exponents as for standard DP in 1 -|- 1 dimensions. 



3 BoundEiry activity and backbone 



Now consider the activity on the absorbing boundary. Direct measurement (by 
box counting) of the fractal dimension of living sites along the time direction 
yields a dimension Di = 0.578±0.005. This is an accordance with the result in 
figure 3 of an independent measurement of the distribution of time intervals 
in which the boundary is dead, i.e., the distribution of first return Pfirst(^) 
of activity. From figure 3 we obtain Pfirst(^) ~ l/t^"''^* with rgrst = 1-577 ± 
0-005, confirming the scahng formula Tfirst = 1 + -Di- We also obtain that 
the distribution of all return Paii(0 ~ l/i"^^" fulfills Tau = 1 — Di. Without 
an absorbing boundary, the dimension of the activity along an arbitrarily 
longitudinal cut of the DP cluster would be = 1 — /9/z/|| = 0.841. Thus, the 
activity on a cut is much less when DP can evolve only on one side of this cut 
than when the activity can move freely back and forth between both sides. 
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Fig. 4. Average lifetime as function of p for values above iii 1 + 1 dimensions. 
The estimate for the slope is ri = 1.000 ± 0.005. 

The connection between DP with a wall and the field-theoretic formulation 
[14] has been worked out in Ref. [19]. Therein it is shown that Pi determines 
the scaling of the density e^^ on the wall. Thus, with the boundary we obtain 
the relation Di — 1 — Pi/v^ — l/i/y in nice agreement with our numerical 
results. 

The backbone is obtained from the infinite cluster by removing all dangling 
ends. Thus the backbone consists of precisely those bonds which would be 
occupied by both the time-directed DP process and its reversed time-directed 
process. It then follows that the backbone density is described by the 
exponent — 2/3. We have numerically confirmed that this backbone be- 
havior is still valid for DP with a wall. Besides, we have measured the backbone 
dimension on the wall with the result Df^ = 1 - l3f^/iy\\ = 0.16 ± 0.01. This 
yields — 1.46 in nice agreement with /3f^ — 2/3i, cf. the result for the 
bulk. 

On the backbone one can identify the so-called red bonds [20] which, if one is 
cut, divide the cluster into two parts. A renormalization group argument [21] 
(see also [11]) yields for the number of red bonds up to scale t the scaling 

NR{t)r^t'/''K (8) 

We have measured the scaling of red bonds for DP with a wall and obtained 
results in complete agreement with Eq. (8). In addition we have measured 
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the scaling of the red bonds along a longitudinal cut for DP with no wall 

and obtained the result N^\t) ~ t-0-04±o.02^ r^j^j^ 

is in accordance with the 

expected behavior N^^(t) ~ t^/'^\\-^ ^-o.o56^ where the extra factor originates 
from the scahng of the width of the cluster. Finally, we have measured the 
scahng of red bonds on the wall for DP with the wall and found the behavior 



4 Results in 2 + 1 dimensions and concluding remcirks 

The above procedure can be applied to higher dimensions. For 2 + 1 dimen- 
sional bond directed percolation on a bcc lattice we have Pc = 0.287338, 
(3 = 0.584, z/|| = 1.295, iy± = 0.734, 6 = 0.451, and x = J^\\/j^± = 0.567 [22,23]. 
We have measured the lifetime distribution, cluster distribution, and average 
lifetime and cluster distributions. Again, we conclude that i/y and i'± are un- 
changed. Furthermore, (3 is needed in order to describe the scahng of the mass 
whereas a new exponent, f3i, is needed to describe the lifetime distribution: 
we find the values j3i = 1.05 ± 0.05, and ri = 0.26 ± 0.02 in accordance with 
Ti = — (3i, cf. Eq. (1). In contrast to the 1 + 1 dimensional case, in 2 + 1 
dimensions (3i does not seem to be (simply) related to the other exponents and 
we conclude that there arc four independent exponents. For 2 + 1 dimensional 
systems with an edge (which introduces yet another exponent) see Ref. [19]. 

One may speculate what the effect of the boundary will be if the cluster is 
initiated some fixed distance x away from the boundary. In that case, scaling 
will be dominated by the DP behavior (1) until the meandering of the cluster 
allows it to reach the absorbing boundary, i.e., until a crossover time ty^ oc 
x'^\\/'^-L_ After this crossover time, the distribution of cluster survival times will 
become steeper, and be determined by (2). Thus the presence of a boundary 
will cause two scaling regimes for the cluster size. One can also consider the 
case where the DP process is initiated in between two absorbing boundaries. 
In such a case the DP process cannot evolve infinitely. After some time, of the 
order of L'^ii^'^^, where L is the distance between the boundaries, the otherwise 
critical DP process dies out. 

In summary, we have investigated the effects of an absorbing boundary on DP 
and discussed the novel scaling behavior originating from such a boundary. 
The scaling can be understood in terms of the three exponents for DP plus 
one new exponent. Our results are in agreement with recent series expansion 
estimates. We also study the scaling at the boundary. We stress that boundary 
effects are not limited to the case where the boundary is absorbing — it could, 
e.g., also be reffecting. The effect that reduces the lifetimes of the clusters and 
the fractal dimension of activity on the boundary is solely due to the fact that 
living branches cannot meander back and forth over the boundary. Finally, the 
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exponent which determines the hfetimc distribution in 1 + 1 dimensions is l/i^y, 
which describes the scahng of the red bonds. In the context of self-organized 
interface growth determined by DP paths there exists a relation which involves 
1 / v\\ and dictates that the corresponding average 'hfetime' should scale as (4) 
with an exponent equal to unity [10]. It would be interesting to explore these 
similarities and to possibly relate the lifetime distribution to the scaling of the 
red bonds. 
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